Viscous cavity damping of a microlever in a simple fluid 
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We consider the problem of oscillation damping in air of a thermally actuated microlever as it is 
gradually approached towards an infinite wall in parallel geometry. As the gap is decreased from 
20 iim down to 400 nm, we observe the increasing damping of the lever Brownian motion in the 
fluid laminar regime. This manifests itself as a linear decrease with distance of the lever quality 
factor accompanied by a dramatic softening of its resonance, and eventually leads to the freezing of 
the CL oscillation. We are able to quantitatively explain this behavior by analytically solving the 
Navier-Stokes equation with perfect slip boundary conditions. Our findings may have implications 
for microfluidics and micro/nano-electromechanical applications. 

PACS numbers: 47.61.Fg, 47.15.Rq, 85.85.+j, 07.79.Lh 
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Micro- and nano-scale mechanical levers are increas- 
ingly used as sensors and actuators in a large variety 
of fundamental studies and applications. Mass detec- 
tion at the zeptogram scale [1], sub-attonewton force 
detection [2] and optical cooling of microlevers [3] are 
among the most spectacular achievements of oscillating 
cantilevers (CLs). These realizations mainly rely upon 
the extraordinary high quality factors (Q) of oscillat- 
ing CLs in vacuum and/or cryogenic temperatures where 
values exceeding 100 000 are attainable. Clearly, main- 
taining such performances in air or in a liquid is a very 
challenging issue as oscillation damping in the surround- 
ing fluid dramatically degrades Q. This has been par- 
tially circumvented by using ultrasmall self-sensing nano- 
electromechanical systems (NEMS, i.e., actuated me- 
chanical devices made from submicron mechanical com- 
ponents facing each other) operating in ambient condi- 
tions of temperature and pressure [4]. 

However, oscillating CLs are also used in viscous envi- 
ronments on many occasions [HI, @, 0, H] • In Atomic Force 
Microscopy (AFM) for example, a resonant CL is used to 
measure surface topography and physico-chemical prop- 
erties of various materials not only in air but also in liq- 
uids for, e.g.^ visualizing dynamic biomolecular pro- 
cesses at video rate [10]. The interaction between an 
AFM CL and a surrounding liquid has been used for 
a distance calibration in a Casimir force measurement 
JjJ and has led very recently to the spectacular demon- 
stration of a repulsive Casimir force [12]. Therefore, the 
need for a quantitative study of the CL behavior in vis- 
cous micro- and nano-scale environments is increasing. In 
this letter, we report such a quantitative study and show, 
down to the submicron scale and in the demanding plane- 
plane geometry, how confinement and boundary condi- 
tions at the solid-fluid interfaces conspire to change the 
coupling to thermal bath and how this can freeze out the 
lever oscillation. 

When a CL beam vibrates in a viscous fluid, the fluid 
offers resistance to the beam displacement If 
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FIG. 1: Scheme of the experimental setup (not to scale). The 
analyzed mechanical system is a standard AFM CL. A cleaved 
optical fiber (bottom) is used as mobile substrate forming a 
cavity with the CL. An etched optical fiber (top) is used for 
the interferometric detection of the CL Brownian motion. The 
main geometrical parameters used in the text are identified, 
the z origin is taken on the flat cleaved-fiber surface. 



the CL is vibrating close to a solid surface, the behav- 
ior of the fluid and, consequently, that of the lever are 
modified by the surface due to confinement. The Navier- 
Stokes (NS) equations give a complete description of the 
fluid behavior taking into account the particular environ- 
ment under analysis. However an analytical solution of 
NS equations is possible only for a restricted number of 
geometries and comparison of theory with experimentally 
relevant configurations is in general a complex matter or 
is even lacking, especially at the deep micron and sub- 
micron scales [l5|, E [l7|, |l8|, 0, H where boundary 
conditions at the fluid-solid interfaces are strongly modi- 
fied [21, 22, [l^,!!^]. In this work, we focus on the dynami- 
cal behavior of a microlever close to a planar rigid surface 
in the air. Provided that adapted boundary conditions 
are used, the NS equations can be solved analytically for 
this plane-plane model geometry that mimics a basic part 
of a MEMS (the counterpart of a NEMS in the micron 
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range) device operating in the air. This, combined with 
the use of the fluctuation-dissipation theorem and of an 
experimental arrangement specially designed to gain ac- 
cess to the intrinsic behavior of the CL, enables us to 
make a quantitative comparison between theory and ex- 
periment in a wide range of cavity lengths down to a few 
hundreds of nanometers. 

Our setup is shown schematically in Fig. [H Its first spec- 
ification is that the CL - a commercial thin silicon AFM 
CL [25I for liquid imaging with dimensions L x w x t = 
107 X 30 X 0.18 jum^ - is actuated by the stochastic ther- 
mal noise only. This induces sub- Angstrom oscillations 
at the CL resonance frequency (cc;o/(27r) :^ 49.5 kHz), 
thereby allowing us to consider the fluid in the cavity in 
the laminar regime. Second, the planar rigid surface fac- 
ing the CL to form a parallel-plate cavity is made of a 
cleaved optical fiber with a diameter of 125 /im that is 
mounted over a three-axis inertial motor so as to be able 
to adjust the cavity gap. This positioning system offers 
a large displacement range (8 mm each axis full range) 
with a good accuracy (40 nm per step). Finally, the CL 
Brownian motion is measured by means of a non invasive 
interferometric detection based on the use of a very thin 
optical fiber facing the CL at a 2 /am distance. This fiber 
has been chemically etched so as to reduce its diameter 
to 5 /im. This corresponds basically to the fiber core di- 
ameter plus a residual amount of the optical cladding for 
better light guidance. The large ratio in excess of 600 
between the areas of the cleaved and detection fibers in- 
sures that only the cleaved one induces air confinement, 
not the etched one, which is used for detection purpose 
only. Therefore, no additional uncontrolled confinement 
and damping are produced by the detection fiber. 
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FIG. 2: The experimental Brownian oscillation power spec- 
trum of the microlever for different cavity gaps. 

An AFM CL vibrating in a viscous fluid may be viewed 
as a driven and damped ID harmonic oscillator whose 
equation of motion reads 



where m, k are the CL effective mass, time- 

dependent position, and stiffness, respectively, 7 is the 
damping factor and F^xt the external (ie., thermal) driv- 
ing force. According to the fluctuation-dissipation the- 
orem, the thermal Brownian motion of the CL at tem- 
perature T is accounted for by a frequency independent 
force power spectrum defined as Sf{^) = 2/cb^7 {^b is 
the Boltzmann constant, 00 the pulsation linked to the 
frequency / = Cc;/(27r)). Starting from Eq. [T] we ob- 
tain the CL displacement power spectrum as Sz{uj) = 
Sf{^)\x{^)\'^ ^ where the CL transfer function x(^) is 
given by: 



m {uq — co'^) — ijco 



(2) 



with ujQ = J ^ . In the limit of small damping, i.e., 

^ <C cjo, Sz{u^) has a resonance at loq. 
In Fig. [21 the experimental Brownian oscillation power 
spectrum 2Sz{u^) is presented as a function of frequency 
/ for different cavity gaps d. It is clearly seen that the res- 
onance peak dramatically broadens and softens to lower 
frequencies with decreasing gap. Within the experimen- 
tal accuracy, we find that the area under the resonance 
curves in Fig. [2] remains constant and equals to the ther- 
mal energy. This shows that the CL damping increases 
with decreasing gap. 

Now, we turn on to a quantitative analysis of the exper- 
iment. The fluid responsible for the CL damping is the 
air confined between the CL and the mobile fiber. The 
dynamic of such an incompressible fluid is described by 
the NS equations 



dv 

m 



- V • Vv 



■ r]V V — Vp, 



(3) 



mz{t) +7i(t) +kz{t) 



(1) 



where v is the fluid velocity, p its density, 77 its dynamical 
viscosity and p the gas pressure. In the laminar regime, 
i.e., in the limit of small Reynolds numbers, Eq. [3] sim- 
plifies to r]V'^v Vp. In order to solve the NS equations, 
one needs to know the specific boundary conditions exist- 
ing at the fluid-solid interfaces (for simplicity we assume 
the cavity plates to be infinitely extended). While for 
macroscopic hydrodynamic applications one usually ac- 
cepts that fluids do not slip against solid walls, this is gen- 
erally not true for microfluidic problems involving MEMS 
or NEMS [21]. A critical parameter in this respect is 
the Knudsen number [2l| Kn = X/d which depends on 
the gas mean free path A. For air at ambient conditions 
A :^ 60 nm which leads here to Kn 0.001 - 0.06. In 
this range of Kn values, it is already known that fluid-slip 
can occur over a solid interface [2l|, [22I, H^] . In particu- 
lar, partial fluid-slip has been recently observed in the 
plane-sphere geometry in air using an AFM in dynamic 
mode [2^. However, none of the previous works inves- 
tigated the regime of Brownian oscillations with typical 
CL amplitudes 5z ~ 0.05 nm much smaller than A {i.e., 
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\/5z ^ 10^). In such a regime, boundary conditions are 
expected to be even more strongly modified [21] com- 
pared with the macroscopic regime although it is not 
yet known how much they are modified. Here, we make 
the hypothesis of perfect slip, historically anticipated by 
Navier [2^, for which friction along the solid interface 
is prohibited, and show that we can obtain a consistent 
quantitative description of our experimental data. This 
hypothesis results in a velocity gradient along the z di- 
rection which leads to a Stokes friction coefficient: 
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2r]A 



(4) 



where A = wL is the cantilever surface. The usual no- 
slip condition at the fluid solid interface would predict 
7 :^ rjwL^ /d^, in clear disagreement with the experiment. 
As a direct consequence, Eq. 4 leads to a much smaller 
decay of the friction force with d than predicted usually. 
Quantitative information on the damping factor is ob- 
tained from the analysis of the CL quality factor. Both 
quantities are linked together by the relation Q 
which becomes for small gaps: 



k 

a;o7 ' 



Q = 7^ jd. (5) 

Eq. [5] predicts a linear dependence of Q with d that can 
be compared with experiment. 
Fig. [3] depicts the quality factor as function of the cav- 
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FIG. 3: The quality factor as a function of the cavity gap. The 
black dots are the experimental points. The insert depicts a 
zoom of the small gap range. In both cases, the red curves 
exhibit the theoretical prediction based on Eq. [5l i.e., the 
prediction of the NS equation in the perfect parallel plate 
geometry. The blue dash-dotted line in the insert shows the 
prediction of the NS model taking into account a residual 
angular misalignment of the cavity as discussed in the text. 

ity gap. Two different regimes can be distinguished. For 
large gaps above 40 /im, the quality factor remains con- 
stant. This is the unconfined fluid regime where no addi- 
tional damping can take place with decreasing gap. For 



smaller gaps however, the quality factor tends to decrease 
with a decreasing gap. We will focus below on the small 
gap regime where the hypothesis of infinite planes is phys- 
ically justified. Since the AFM CL has a surface 10 times 
smaller than the substrate fiber, the gap limit for the hy- 
pothesis of infinite planes to remain physically sound can 
be estimated by taking the apparent CL surface as refer- 
ence, i.e. d^Y^ ~ \h\im. 

As shown in the inset of Fig. [3] the experimental results 
and the theoretical prediction of Eq.[5l with no adjustable 
parameter [2^, coincide to within 5% for gaps larger than 
5 livn but for smaller separations, the agreement worsens 
to reach 100% at the smallest gap, 400 nm. We interpret 
this difference with a residual small angular misalign- 
ment of the two facing parallel plates. For small misalign- 
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FIG. 4: The resonance frequency as function of the cavity 
gap in the small gap regime. Like in FigU the red dashed, 
respectively blue dash-dotted curves is the prediction of the 
NS model for the perfectly aligned, respectively slightly mis- 
aligned, cavity. 

ment, the problem can be treated within an approxima- 
tion similar to the Proximity Force Approximation (PEA) 
used, for instance, in the Casimir force formulation in the 
sphere-plane geometry [28]. In this approximation, the 
corrected damping factor becomes: 
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dx • dy 



do -\- X tan a -\- y tan /3 



(6) 



where do is the shortest distance from the inclined CL to 
the substrate, a and P are the lateral tilt angles of the CL 
with the mobile surface in the x and y directions, respec- 
tively. The angle values that permit to reproduce the evo- 
lution of the disagreement between experiment and the- 
ory, as presented in Fig.[3l are a ^ /3 ^ 10 mrad. Consid- 
ering these misalignment angles, the good agreement be- 
tween theory and experiment can now be extended down 
to the smallest gap range that we have measured as can 
be seen in Fig.[3l Over the entire range 400 nm — 15 /im, 
the remaining disagreement is ^ 5% only [29j . 
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We now discuss the frequency softening of the CL os- 
cillation, the other salient experimental fact revealed by 
Fig. [21 In the limit of large damping, i.e. the approx- 
imation ^ <C no longer holds, the power spectrum 
of Eq. [2] has a down-shifted resonance pulsation uj' given 
by; 

Fig. m shows that the resonance frequency shift {uj' — 
cJo)/(27r) can be extremely large. One step further below 
400nm, a complete freezing of the CL oscillation would 
have been observed. This was precluded by the resid- 
ual angular misalignment discussed above. Taking into 
account for data analysis of the misalignment obtained 
from Fig. O we can quantitatively model the measure- 
ments in Fig. m without any adjustable parameter what- 
soever. Therefore, beside the CL width that governs the 
Q factor behavior at large scale (Fig.[3|), another shorter 
characteristic length is here emphasized in the submicron 

range, i.e. c^^rit ~ "^^^ (around 500 nm in our case), 
which is the gap width cancelling the resonance frequency 
in Eq. [71 This characteristic length is determined by the 
CL dynamics and the fluid viscosity. 
In conclusion, we have presented high sensitivity mea- 
surements of the damping of a thermally driven CL in a 
simple fluid confined in a microcavity formed by this CL 
facing an infinite wall. As the cavity length decreases, the 
fluid confinement induces a dramatic damping of the CL 
Brownian motion which can lead to its complete freezing 
at small gaps. A consequence of our work is that micro- 
or nano-oscillators can either present high Q factors or 
be overdamped systems depending on their actual geom- 
etry, resonance frequency, oscillator substrate gap and, of 
course, ambient viscosity. These findings may impact the 
design of modern NEMS and microfluidic devices since 
the 1/d dependence strongly reduces dissipation even for 
separations d as large as thousands of mean free path A 
(see Fig. 3). This 1/d behavior can be furthermore de- 
scribed by solving the Navier- Stokes equation with per- 
fect solid-fluid slip boundary conditions. The agreement 
between experiment and our model is found over a broad 
range of cavity lengths, including the submicron range. 
Interesting extensions of the present work include the 
study of parameters affecting boundary conditions, such 
as external actuation of CLs (to obtain large oscillation 
amplitude) [13] jUanostructuration [4] , and surface chem- 
ical properties [8]. 
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